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The radial alignment of dark matter subhalos: from 
simulations to observations 



oo ■ 

o ■ 

O 



I 

o 



> 

cn 

00 

in 
o 
oo 
O 



X 



Alexander Knebe^, Hideki Yahagi^'^'^'^'^, Hiroyuki Kase^'^, Geraint Lewis^, 
Brad K. Gibson^'^ 

^ Astrophysikalisches Institut Potsdam, An der Sternwarte 16, Germany 

^Department of Astronomy, University of Tokyo, Tokyo 113-0033, Japan 

^Division of Theoretical Astronomy, National Astronomical Observatory, Japan 

'^Research Fellow of the Japan Society for the Promotion of Science 

^ Centre for Astrophysics, University of Central Lancashire, Preston PRl 2HE, UK 

^ Research Institute for Information Technology, University of Kyushu, Fukuoka 812-8581, Japan 

"^present address: Prometech Inc., 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 

^School of Physics, University of Sydney, Sydney NSW 2006, Australia 



Submitted Version 



ABSTRACT 

We explore the radial alignment of subhalos in 2-dimensional projections of cosmolog- 
ical simulations. While most other recent studies focussed on quantifying the signal 
utilizing the full 3-dimensional spatial information any comparison to observational 
data has to be done in projection along random lines-of-sight. We have a suite of well 
resolved host dark matter halos at our disposal ranging from 6 x 10^^/i~^Mq down 
to 6 X lO^^/i~^M0. For these host systems we do observe that the major axis of the 
projected 2D mass distribution of subhalos aligns with its (projected) distance vector 
to the host's centre. The signal is actually stronger than the observed alignment. How- 
ever, when considering only the innermost 10-20% of the subhalo's particles for the 
2D shape measurement we recover the observed correlation. We further acknowledge 
that this signal is independent of subhalo mass. 

Key words: galaxies: evolution - galaxies: halos - cosmology: theory ~ cosmology: 
dark matter - methods: A''-body simulations 



1 INTRODUCTION 



The concordance of a multitude of recent cosmological stud- 
ies has demonstrated that we appear to live in a spa- 
tially flat , A-dominated cold dark matter (ACDM) uni- 
verse (cf. ISpergel et aLll200'i1 ). During the past decade sim- 
ulation codes and computer hardware have advanced to 
such a stage where it has been possible to resolve in de- 
tail the formation of dark matter halos and their subhalo 
popu lations in a cosmological context (e.g. iKlvpin et al.l 
1 19991 ). These results, coupled with the simultaneous increase 
in observational d ata (e.g . 2 deg ree Field galaxy redshift 
survey (2dF GRS). ICoUesj (12003) ; Sloa n Digital Sky Sur- 
vey fSDSS). lAdelman-McCarthvl l|2007^ ). have opened up 
a whole new window on the concordance cosmogony in 
the field that has become known as "near-field cosmology" 
iFreeman fc B land-Hawthorn 2002). 

One particular property of the satellite population that 
has caught the attention of simulators recently is the radial 
alignment of their primary axes of subhalos with respect to 
the distance vector to the host's centre. The first evidence 



for this effect was reported for the Coma cluster, where it 
was observed that the projected major axes of cluster mem- 
bers pr eferentially align with the direction to th e cluster 
centre iHawlev fc PeeblesI 1 19751 : iThompsonI 1 19761 ). Such a 
correlation between satellite elongation and radius vector 
has further been confirmed by statistical analysis of the 
SDSS data jPereira fc Kuhnll2005l: lAgust sson fc BrainerdI 



l2006l : iFaltenbacher et all 120071 : IWang et al.ll2007l ). The 

dial alignment of subhalo shapes towards the centre of their 
host has also been measu red for the subhalo population in 
cosm o logical simulations llKnebe et al.l l2008l: iKuhlen et al.l 



I2OO7I : IFaltenbacher et al.l I2OO7I : iPereira et all |2007| ). We 

though note that all these authors used the 3-dimensional 
spatial information available to them. However, any fair 
comparison to the signal found observationally requires the 
projection of the simulation data into an observer's plane, 
i.e. averaging ov er a substantial number of 2D restrictions 
of the data (cf. IFaltenbacher et alll2007l ). This is the pri- 
mary motivation for the present study: Can we still find a 
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signal of radial alignment when limiting the simulation data 
to 2D? 

In this Letter we provide evidence that the radial align- 
ment of subhalos in cosmological simulations - when pro- 
jected into two dimensions - is substantially stronger than 
found in observations. We though recover the observed cor- 
relation strength when restricting the shape measurement 
to the innermost 10-20% of the sub halos' partic l es; th is re- 
sult supports the findings of iFaltenbacher et alj (|2007l ) who 
also used the inner regions of their satellite halos as a proxy 
for the orientation of a hypothetical galaxy and found good 
agreement with the measurements from SDSS data. We fur- 
ther show that the signal does not depend on the mass of 
the actual subhalo. 



2 METHOD 
2.1 Simulations 

The cosmological A'^-body simulation employed in this study 
was g enerated as part of the N umerical Galaxy Catalog pro- 
gram l|Nagashima et al]|2005l). Utilizing a parallel version of 
an es tablished A''-body code l|Yahagi fc Yoshiill200ll : lYahagil 
I2OO5I ). the algorithm uses an adap tive mesh refinement 
technique to main t ain accuracy (cf. Kravtsov et al.' ' 19971: 
iKnebe et all I2OOII : lYahaei fc Yoshiil I2OOII : l^yssior 20o3). 
The simulation considered a region of the universe with a 
boxsize of 70 Mpc (comoving), with 512"^ particles, re- 
sulting in a particle mass of 3.04 X lO*/i"^M0. The adopted 
cosmological parameters represent a concordance ACDM 
model with Qrn = 0.3, = 0.7, h = 0.7, and as = 0.9 (cf. 
ISpergel et a"Lll2003l ). The simulation was started at z = 41 
and analysed at z = 0. 



Table 1. Number of subhalos in a certain mass range. 



2.2 (Sub-)Halo Identification 

In extracting groups and subhalos, the standard friends - 
of- friends (FoF) method was employed (jPavis et al.lll985h . 
with the ratio of the linking length to the mean sepa- 
ration of particles b = 0.2, and cluster and group scale 
host halos were choosen to have masses is in the range, 
6.45 X IO^^/i-^Mq^ Mhoat ^ 5.95 x IO^/i-^Mq. Once these 
host halos were identified , a hierarchical fin der based on 
the FoF (HFoF) method ijKlvpin et all Il999l ') used to iso- 
late subhalos, employing s successively decreasing linking 
length. Furthermore, the eva porative method was used to 
discard unbound subhalos (jPfitzer et al.|[l997l ). and all un- 
bound particles were removed from the subhalos. At this 
stage, we calculated potential energy of particles in subha- 
los iteratively, defining a minimum subhalo as containing 100 
particles. Q While this is very time consuming procedure, the 
process was accelerated by employing a special purpose com- 



puter for self-gravitating system, GRAPE-6 (|Makino et al 
I2OO3I) : details of the subhalo finding are given in lKase et al 
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We ended up with 40 host halos with a combined num- 
ber of 2188 subhalos in excess of 100 particles. In the fol- 
lowing analysis we stack the information from all these hosts 
and present the results in one single plot. We are confident 
not to obscure any signal as it has recently been shown 
that the radial alignment is independent on the host mass 



that the radial alig r 
l|Knebe et al.ll2008l i. 



When investigating the dependence of radial alignment 
on the mass of the subhalos we split the subhalos into three 
mass bins, i.e. [lO", 10^^]/i-^Mq, [10^^ 10^^]/i-^MQ , and 
[10^^, 10^'^]h~^MQ. The total number of subhalos in each of 
these bins (as well as the number of subhalos with M < 
lO"/i"^M0 is given in Table[T] We further like to note that 
the most massive subhalo weights M^^" = 2.9 x lO^/i'^M© 
and belongs to a A/host ~ 1-6 x 10^'^h~^MQ host halo. 



3 RESULTS 

Our view of clusters of galaxies upon the sky is a projec- 
tion of the true, three dimensional distribution and orienta- 
tion of halos and their associated galaxies. In this study, we 
therefore consider projected views of our simulated halos, 
adopting 19 projections in both azimuthal angle and cosine 
of latitudinal angle, respectively. The results presented sub- 
sequently in this study average over these differing orienta- 
tions. 



3.1 Subhalo Shapes in 2D 

The principal axes of the projected two-dimensional den- 
sity distributions of all our subhalos were ca lculated us- 
ing the reduced moment of inertia tensor (|Kat j Il99ll : 
iDubinski fc Carlberd[l99ll i^. 
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where xi^i is i-th component of projected (2-dimensional) 
position of the Z-th constituent particle, and is the num- 
ber of particles in the subhalo. Solving for the eigensystem 
of this tensor identifies the principle axes of the distribution. 
We then define the sphericity 
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^ We note that this is a conservative criterion as other studies 
indicate that halo catalogues are complete for objects containing 
in ex cess of ^50 particles (e.g. iKravtsov et al..2004 : .Allgood et al.l 
I2OO6I) . 



We like to note in passing that it makes little difference using 

either the reduced or the "standard" moment of inertia tensor 
/-2 



that does not include the 1/C; term. 
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Figure 1. Distribution of the 2D sphericity s = a^jax. Solid 
line indicates the distribution of axis ratio using all subha- 
los. The different lines represents restrictions of the subhalo 
masses under consideration to the intervalls [10^^, 10^^]h~^MQ, 
[10^2, lol3]/i-iM0, and [lO", IQi^J/i-iMg. Error bars indicate 
the square root of the sum of the Poisson error and the projection 
variation. 



of a (projected) subhalo as the ratio of the minor to ma- 
jor axis. The procedure in determining these principle axes 
has to be done iteratively, starting with the assumption of 
sphericity s = 1. The eigensystem is solved for s until the 
value converges. 



3.1.1 Mass Dependence 

In Fig[l] we examine the variation of the (normalized) 
probability distribution P{s — 02/01) of sphericity with 
subhalo mass. We observe a (marginal) trend for s to 
decrease with incre asing subhalo mass - as reported by 
lAllgQod et al.l ()2006l ) for field halos. We also note that the 
projected spheric ities are o f similar value to the 3D spheric- 
ities as found in .Knebe et al.. (2008. ') making a study of ra- 
dial alignment feasible. We need to stress that the num- 
ber statistics for the subhalo mass bin [10^'^, 10^*]/i~^Mq 
is rather small (cf. Table [1]); there are a mere 23 sub- 
halos in the respective mass range and hence the large 
error bars. As these objects are in fact rather massive 
and comparable to the actual host halos, we are not sur- 
prised to find sphericities closer to those o f isolated/field 
halos ((s) » 0.66, cf. iFrenk et al.1 Il988l: IWarren et all 



19921: ix asun fc Evrar dI l2005l: 'Bailin fc Stcinmct: 



I I I 

z 2005; 



Allgood et al..,2006i : iMaccio et al.,,2007. : .Bett et al...2007. ') 



3.1.2 Radial Dependence 

As with many cosmological A''-body studies, we are con- 
sidering the orientations of dark matter (sub-)halos and 
are comparing them with the distributions of baryons ob- 
served in galaxy clusters. Such a comparison can be signif- 
icantly biased if the observed orientation of galaxies, which 
occupy the deepest part of the subhalo potential, are not 
aligned with the global o rientation of the dark matter (cf. 
iBaihn fc Steinmetz|[2005l ). Hence we also present the varia- 
tion of subhalo sphericity (as well as the radial alignment in 
Section [3. 2. 2|) with distance to its centre. 



Figure 2. Distribution of the 2D sphericity s = 02 /ai. The dif- 
ferent lines represent shape measures at different mass thresholds. 
Error bars indicate the projection variation. 



To select the inner parts of each subhalos, (^i was calcu- 
lated assuming s = 0. Particles were then sorted by Q and 
the inner subhalo was selected as a numerical fraction of 
the subhalo's mass. Once this distribution was defined, the 
moment of inertia tensor was recalculated, and the eigensys- 
tem solved iteratively to determine the inner value of s (as 
outlined in Section [3. ip . 

FiglH shows that the inner regions of subhalos appear 
to be significantly more elliptical than their overall matter 
distribution. This is in agreement with the findings of others 
who also reported that dark matter hal os tend to increase 
their asphcricity towards the centre (e.g. Ijing fc Sutoll2002l : 
Allgood ct al. 2006; Hayashi et al. 2003). We will return to 
this finding in the following Section when investigating the 
radial alignment as a function of distance to the centre: Is the 
signal actually enhanced due to the increase in sphericity? 
Or does the correlation weaken despite this result? 

We tough like to note that lAllgood et all l|2006l ') also 
cautioned that the sphericity may be overestimated if too 
few particles are used. In order to check the credibility 
of the trend seen in FigH] we recalculated the distribu- 
tions restricting the haloes to have masses in the range 
M e [10^^, lO^^]/i"^M0. We still observe the same shift in 
the peak of the distributions; however, the extended tail to- 
wards lower s-values for the inner 10-20% has vanished, i.e. 
there are hardly any objects with s < 0.4. 



3.2 Radial Alignment in 2D 

The principal aim of this study is to investigate whether or 
not there is a dependence of the radial alignment of subhalos 
(i.e. the alignment of their major axis with respect to the 
centre of the host) on the mass of their host halo when ob- 
serving the mass distributions in projection. Recent observa- 
tional evidence suggests that the major axis (in projection) 
of satellite ga laxies tend to "point towards the centre of their 



host" (e.g. 



[ Pereira fc KuhnI 20051 : lAgustsson fc BrainerdI 

I2OO6I ; iFaltenbacher et al.l 120071 ). It is therefore natural to 
ask whether or not subhalos in cosmological simulations 



tigated this subject (iKnebe et al.l 20081: 


Kuhlen et al.ll2007l: 


iFaltenbacher et al.l l2007l: Pereira et al.l 


2007f) all of which 



used the full 3D spatial information available to them; 
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Figure 3. Correlation of angle between major axes of subhalos 
and position of the subhalos from the center of their host. We use 
the same subhalo mass bins as in Fig[T] Error bars (marginally 
shifted for clarity) are again the square root of the Poisson error 
and the projection variation. 



onlv lFaltenbacher et al.l (|2007t ) also presented a 2D analysis. 
However, a fair comparison to observational findings requires 
a projection of the simulation data into the observer's plane. 

To measure the radial alignment of our (projected) sub- 
halos, we use the eigenvector Ea^ which corresponds to the 
direction of the major axis oi of the subhalo. We quantify 
the radial alignment of subhalos as the angle between the 
major axis E^j of each subhalo and the radius vector of the 
subhalo in the reference frame of the host: 

sub 



|E„ 



(4) 



3.2.1 Mass Dependence 

In Fig|3]we show the normalized distribution of the angle (j) 
as defined via Eq.(|3]). While the black line shows the signal 
for all subhalos under consideration, the different lines are 
representative for various mass bins. We find a positive ra- 
dial alignment sign al different from is o tropy, in agreement 
with other studies (Kncbo ct al. 2008; K uhlen et alll2007l : 
iFaltenbacher et "al., 2007: Pereira et al. 200?! ) - even when 
restricting the analysis to 2D. We further conclude from 
this plot that no particular mass range is responsible for 
the signal. We also note that the strength of the signal is 
substatially stronger t han the one seen in observations (cf. 
IPereira fc Kuhnll2005r ). We will though see in the following 
Section that this sensitively depends on the point where we 
measure the actual (projected) shape of the subhalo. 



3.2.2 Radial Dependence 

As with many cosmological A''-body studies, we are consider- 
ing the orientations of dark matter halos and are comparing 
them with the distributions of baryons observed in galaxy 
clusters. Such a comparison can be significantly biased if the 
observed orientation of galaxies, which occupy the deepest 
part of the subhalo potential, are not aligned with the global 
orientation of the dark matter. We examined this effect by 
studying the dependence of the (projected) radial alignment 
signal on the point where we measure the subhalos' shape 
as outlined in Section [3. 1.21 
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Figure 4. Same as Fig[3l but each line shows the correlation using 
major axes derived from a fraction of particles. The different lines 
represent shape measures at different mass thresholds. Error bars 
indicate the projection variation. 



The result of selecting the inner particles is quite strik- 
ing; Figl4]presents the orientation correlation using different 
percentage of inner particles in subhalos. We observe a drop 
in the signal's strength when restricting to the central parts; 
considering only the inner 10-20% of particles we recove r 
the obser ved correlation not e d by Pereira fc KuhnI (|2005|): 
follow ing IStruble fc Peebles! (|l985l ) and iPereira fc KuhnI 
(|2005l ). respectively, we quantify our results in terms of an 
alignment parameter 5 = '^.{4>i/N) — 45°, where (j>i are the 
indivudal subhalo orientations and A'^ the number of subha- 
los. When using all particles in a subhalo we derive a value 
of 5 ~ —11°. However, when restricting to the innermost 
20% (10%) we obtain 6 ^ -2° (S » -1°) in a greement with 
the value reported by IPereira fc KuhnI l|2005l ). The weaken- 
ing of the signal could be due to poorly defined principle 
axes if the inner distributions of particles were signficantly 
rounder than the overall dark matter distribution. However, 
we already noted that the distribution of sphericity moves 
towards smaller values in the inner parts of the subhalos (cf. 
FiglH)- In that regards caution must also be urged as, for 
the smallest halos containing of order 100 particles (cf. Sec- 
tion (521), considering the inner 10% of particles means that 
we are measuring principle axes using of order 10 particles, 
and as such are very sensitive to Poisson noise. In order to 
avoid such small numbers of particles and to verify the credi- 
bility of the weakening of the signal, we repeated the analysis 
restricting the subhalos to those whose mass is in the range 
lO^^/i"^M0^ Msnb < lO^^/i"^M0, that is the number of 
constituent particles is in the range, 3289 ^ Np < 32897. 
We recover the same orientation correlation, with the lower 
the fraction of inner particles, the lower the observed sig- 
nature of orientation correlation, again bringing into agree- 
ment with the observed correlation for the inner 10-20% of 
particles. 



4 SUMMARY AND CONCLUSIONS 

We have examined whether or not the radial alignment of 
the major axes of subhalos with respect to the centre of 
their host dark matter halo is still present when projecting 
the data from cosmological simulations into the observer's 
plane. Our results draw upon a sample of dark matter host 
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halos spanning group- to cluster-mass scales (6xlO^^/i"^M0 
to 6 X lO"/i"^M0). 

While observations of clusters of galaxies have revealed 
a mild radial alignment of subhalos with r espect to their host 
l|Pereira fc Kuhnl2005l : IWang et al.ll200'i1 ) the signal found in 
our data is substantially stronger. However, when restrict- 
ing the shape determination to using only the innermost 
10-20% of the subhalo's particles we recover the observed 
co rrelation strength, in agree ment with the results reported 
bv iFaltenbacher et al.l l|2007h . We though note that the sig- 
nal does not weaken because the central parts of subhalos 
are rounder and hence radial alignment itself becomes ill- 
defined. On the contrary, in agreement with previous find- 
ings we observe the trend for subhalos t o become more as- 
pherical when approaching the centre (cf. Ijing fc Sutol200^ : 
lAllgood et aLll2006l : iHavashi et al.ll2007l ). Therefore, while 
the strength of the radial alignment decreases even though 
the asphericity increass towards the centre there has to be 
a misalignment between the inner regions of the subhalo 
and its outer structure. However, it is known that the major 
axes of the central parts of field halos are aligned with those 
of the whole object (cf. Jing & Suto 2002;). One possible 
solution to this puzzling observation is tightly coupled to 
the overall expla nation of the radial alignment as a dyn am- 
ical/tidal effect (|Kuhlen et al.ll2007l : [Pereira et al.ll2007^ : it 
can be surmised that subhalos which have passed through 
their apocenter must experience a strong tidal field and the 
outer - more loosely bound - regions of a subhalo will be 
subject to more distortion than the inner parts. This results 
in a weaker correlation of the protected inner parts parts of 
the sub-halo and the host. 

We may even be as bold as to reverse the argumentation 
and use the (projected) radial alignment to infer the point 
where the galaxy is expected to lie within the dark matter 
subhalo: The point where the strength of the signal matches 
the observed correlation - that is based upon the shape of 
the galaxy - coincides with 10-20% of the total subhalo mass. 
As suming a dens i ty pro file of the functional form proposed 
bv lNavarro et all l| 19971 . NFW) and a fiducial concentration 
of about c=10-15 for this profile, 15% of A/vir is reached 
at about 9-11% of the viri al radius |f| This is i n agree ment 
with the results found by iBailin fc Steinmet j l|2005f ) who 
found that this marks a point where the influence of the 
baryons becomes important and there appears to be a lack 
of connection between the inner and outer regions deflned 
by this point, respectively. 
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